Dynamics of Morphology-Dependent Resonances by Openness in Dielectric Disk for 

TE polarization 
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We have studied the parametric-evolution of morphology-dependent resonances according to the 
change of openness in a two-dimensional dielectric microdisk for TE polarization. For the first 
time, we report that the dynamics exhibits avoided resonance crossings between inner and outer 
resonances even though the corresponding billiard is integrable. Due to these recondite avoidances, 
inner and outer resonances can be exchanged and Q-factor of inner resonances is strongly affected. 
We analyze the diverse phenomena arisen from the dynamics including the avoided crossings. 

PACS numbers: 42.25.-p, 42.55.Sa, 05.45.Mt 



I. INTRODUCTION 

Morphology-dependent resonances (MDRs), which is 
defined as the resonances found in cylindrical, spherical 
and ellipsoidal shaped optical cavities, are a much studied 
topic in the field of optics due to their practical impor- 
tance in many applications [1| as well as their heuristic 
value. More than a decade ago, Johnson is the first to 
approach MDR in a spherical cavity [2| as an analogy 
of quantum mechanical resonances, which is familiar in 
nuclear, atomic, and molecular scattering theory. By the 
same approach, Nockel have studied MDR analytically 
in a cylindrical cavity [3| as an unperturbed system in 
order to analyze a slightly deformed one. 

Recently, study of the MDR is resurfaced as a hot is- 
sue in connection with the application of deformed mi- 
crocavity lasers because of its high potentiality of ap- 
plication for photonics and optoelectronics. Bogomolny 
et al. have reported the existence of a different set of 
resonances, which can be classified as shape {outer) res- 
onances [ORs] with high leakage along with Feshbach 
[inner) resonances (IRs) in a two-dimensional dielec- 
tric circular microdisk for transverse magnetic (TM) and 
transverse electric (TE) polarizations [J, Q. The simi- 
lar topic was treated in a numerical method for quantum 
billiard [y|. The evolutions of these two classes of res- 
onances to the so-called small opening limit where the 
refractive index n of the disk goes to infinity have been 
explicitly studied [3| and the applicability of effective po- 
tential analogy for the ORs has been checked Q. One 
intriguing fact of the MDR dynamics from last two stud- 
ies is that, for TE polarization with magnetic field H 
perpendicular to the cavity plane, both IRs and ORs are 
found in the same region of complex wavenumber plane 
where the regime of n is practically relevant. 

In two-dimensional chaotic quantum billiards, it is well 
known that the level dynamics resembles dynamics of 
one-dimensional particles with repulsive interaction as a 
system parameter varies Q. It is called avoided level 
crossing. Similarly, avoided resonance crossing (ARC) 
phenomenon in open Hamiltonian systems is generally 
categorized into two classes [l3|. In the case of a 
strong interaction, shortly strong ARC, the imaginary 



parts of complex-valued energies are crossing and the 
real parts are avoided. The identities of the two modes 
are exchanged via the ARC and their wave patterns are 
strongly superposed during the ARC. For a weak inter- 
action, shortly weak ARC, the real parts are crossing and 
the imaginary parts are avoided while their identities are 
retained. 



In two-dimensional chaotic dielectric cavities, there 
have been undertaken many theoretical |lll4l6| and ex- 
perimental |17Hig | researches related with the ARC. Due 
to the open property, the resonance dynamics in a di- 
electric cavity is much complicated and shows diverse as- 
pects compared to that of the closed billiard. Recently, 
it has been represented that, when a deformation pa- 
rameter is varied in the rectangular or elliptic cavity, of 
which the corresponding billiard is integrable, ARC can 
occur by the effect of openness |2C)I.|21|.^ Another recent 
report have claimed that, even on a slight deformation 
from the disk shape, a symmetry breaking of escaping 
rays occurs, which gives rise to dramatic changes in emis- 
sion patterns depending on the characteristics of defor- 
mations |22| , and emission directions and resonance pat- 
terns are deeply influenced by the refractive index of the 
cavity, i.e., openness of the system, as well as the under- 
lying classical manifold structures [2J| . These imply that 
analyzing the emission directions and mode patterns of 
slightly deformed cavities requires understanding of the 
MDR dynamics for openness. 



In this paper, we uncover the mechanism behind the 
MDR dynamics in a circular disk for TE polarization as 
the openness parameter Q, = 1/n is varied and observe 
the impacts of pure openness on the dynamics. The per- 
fect rotational symmetry provides the analysis possible. 



^ In the pure one-dimensional system like a double-well potential, 
integrability does not forbid avoided crossing between two states. 
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FIG. 1: (Color online) ARCs between IR and low-leaky OR. 
Solid, dotted, and dashed lines are i()(25,8), io(25,9), and 
oo(25, 13), respectively. 
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FIG. 2: (Color online) Metamorphosis of two IRs and a 
low-leaky OR. Top, middle, and bottom rows are io(25,9), 
Oo(25, 13), and io(25,8), respectively. Five pictures in each 
row are at fi = 0.435, 0.448, 0.465, 0.481, and 0.500. 



II. AVOIDED RESONANCE CROSSINGS IN 
DIELECTRIC DISK FOR TE POLARIZATION 

The generic two-dimensional dielectric cavities have 
mainly two types of changeable system parameters, e.g., 
degree of geometrical deformation and openness. In such 
cavities, IRs can generally interact with each other. Their 
resonance dynamics contains a complicated mixture of 
the aspect by deformation and the effect of openness. On 
the other hand, the circular disk does not have any defor- 
mation parameter and the only relevant parameter is the 
openness. Moreover, its rotational symmetry makes that 
the angular component of the wavef unction is trivial. As 
a result, it renders the system effectively one-dimensional 
and then leads more clear and basic study for the pure 
openness without the aspect by deformation. 

The Helmholtz equation of dielectric cavities is given 
by V^V + (k/i^) -0 = 0, where ip and k are the wave- 
function and the wavenumber, respectively. In the disk 
with the radius i? = 1, the complex- valued wavenum- 
bers kr for resonances are obtained from the boundary 
matching conditions for TM and TE polarizations un- 
der the purely outgoing wave condition [24, 25]. In the 
opening regime < 17 < 1, the IRs, i(m,l), are labelled 
by angular momentum quantum number m and modal 
ordering number I which corresponds to the number of 
radial intensity peaks inside the cavity. Similarly, the 
ORs can be also represented by o{m,l), where I is as- 
signed in accordance with the increasing order of Re[A:r]. 
For a given m, the maximum numbers of ORs in TM and 
TE polarizations are [^j and [^^^^] , respectively (the 
notation [] is integer value). Notice that, in TE case, 
some ORs have relatively low leakage and we have espe- 
cially called them low-leaky outer resonances (low-leaky 
ORs) Q . Since there is only one low- leaky OR having the 
largest Re[A:r-] among the ORs for each m(> 0), low-leaky 
ORs can be denoted by o(m, [^^^])- Additionally, the 
notations \o{m,l) and Oo(?Ti, /) are used for IR and OR 
in the small opening limit, respectively, where they can 



be unambiguously identified 0, |g|. The subscript zero 
means fi — ?► 0. 

When Oo(to, [^^^^^]) is evolved from the small open- 
ing limit toward the large opening regime, both weak 
ARC and strong ARC are unexpectedly observed in the 
dielectric disk. Figure [T] shows the two types of ARCs be- 
tween IR and low-leaky OR with m = 25. At SI ~ 0.448, 
Oo(25,13) interacts with io(25,9) via weak ARC. When 
the openness increases, Oo(25, 13) interacts with io(25, 8) 
via strong ARC at $7 ~ 0.481. The radial intensity pat- 
terns at five points on each line in Fig. [T] are shown 
in Fig. [21 Due to the superposition near the point of 
the weak ARC &i Q. ~ 0.448, the leakage of o(25, 13) 
(Fig. Etg)) decreases while that of i(25, 9) (Fig. E^b)) 
increases. However their identities are retained after 
the ARC as shown in Figs. [SJc) and (h). Likewise at 
the weak ARC, o(25, 13) and i(25, 8) also undergo the 
change of leakage and their spatial patterns are super- 
posed near the point of strong ARC as shown in Figs. 
[21[li) and (m). Especially, just at the point of the strong 
ARC, two resonances are strongly superposed and then 
their patterns are almost identical as shown in Figs. 
[5Ji) and (n). After the strong ARC, the identities of 
i(25, 8) and o(25, 13) are exchanged as shown in Figs. [51Jj) 
and (o). Consequently, Oo(25,13) becomes i(25,8) and 
io(25,8) becomes o(25, 13). On a further increase of Vt 
over 0.51, io(25, 8), which has become o(25, 13), interacts 
with io(25,7) through another strong ARC. Eventually 
io(?7i, I) can experience strong ARCs up to two times over 
the whole regime of Vt. io{?7i,l) changes to 0(771, [^^^^^]) 
at the first strong ARC and io{m, I) to i{m, I — 1) at the 
following second strong ARC: 



io{m,l) "' 



strong ARC 



o(to. 



m + 1 



strong ARC 






1(771, I — 1). 



On the other hand, io(25,9) and io(25,8) do not in- 
teract each other, even if they have the same symmetry 



and their imaginary values of kr cross at fi ~ 0.468. It 
has been known that two-dimensional optical dielectric 
cavities are typically non-integrable systems due to their 
open property and the openness breaks the integrability 
of the systems. Consequently, the avoided crossings be- 
tween the IRs occur in the optical cavities, even if their 
corresponding billiards are integrable. It was explicitly 
shown in elliptic and rectangular cavity [2CJ,[21|. But, un- 
like the cavities which correspond to integrable billiards 
such as rectangular and elliptic cavity, the avoided cross- 
ings between IRs intriguingly do not occur in the case 
of circular disk cavity even though the system is open. 
This fact implicitly indicates that the IRs remain with 
the properties of the normal modes in the corresponding 
billiard which is integrable in this case. 



III. TOY MODEL 
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FIG. 3: (Color online) Analytical modeling of (a) weak ARC 
and (b) strong ARC. Solid and dashed lines are IR and OR, 
respectively. 



the zeros of Bessel function 



Jm.l 



and the zeros of 



Hankel function derivatives (7|, /i^ ;, respectively. Then, 
i?i(ri) and Eo{^) for TE case can be simply expressed as 
follows: 



Before we discuss the numerical results, we will con- 
sider the analytical approach through a toy model. The 
ARCs can be displayed with effective Hamiltonian com- 
posed of a non-Hermitian 2x2 matrix and a coupling 
term as follows: 



H,sin) 



Ei{n) 
Eo{n) 



rin 
(,n 



(1) 



where Ei(fl) and Eo{^) are the complex energy fc^ of IR 
and OR, respectively. The off-diagonal element in the 
coupling term, 77 G C (^ G C), describes the coupling 
effect from OR (IR) to IR (OR) and it is assumed that 
the coupling effects are proportional to fi due to the fact 
that the coupling must becomes zero in the small opening 
limit. Equation llj leads the eigenvalues 



i^,(o) = Ml)±Ml)± 



iEm^3^mr+,^n^, 



(2) 

To simply illustrate the ARCs, we substitute Ei{il) 
and Eo{^) with the intuitive relations of the MDR dy- 
namics. Firstly, to obtain the basic MDR dynamics for 
TE case except for the effects of mode interaction and 
Brewster angle by openness, we borrow the dynamics for 
TM polarization, in which the interactions between the 
two types of resonances are infinitesimal in the experi- 
mentally feasible opening regime due to the absence of 
low-leaky ORs. From this, Re[fcr] and Im[fcr] of an IR 
become almost linear to fl and to 7, respectively [3|-(a|- 
Here, 7 given by 



7 = In — 

' 2 i-n 



(3) 



is the lower bound of Im[fcr] of IRs in TM case. Also an 
OR in TM case maintains almost an constant complex 
valued k^ for the change oi fl [J, [a, Q- Secondly, we 
consider that the destinations of nkr of IRs and kr of 
ORs in the small opening limit for TE case, correspond to 



Eo{n) 



Ku 



2 



(4) 



Even though we borrowed the MDR dynamics for TM 
case, the dynamics of y/EjJT) and yj Eo{^) was well 
matched to the numerically obtained resonance tracing 
behaviors [a| for TE case in the opening regime where 
the interaction is very small. Figures [3l^a) and EJb) 
show weak ARC of io(25, 9) and strong ARC of io(25, 8) 
with Oo(25, 13), respectively, when the coupling constant 
77^ = 12737.8. The value of coupling constant was prop- 
erly chosen as compared with the numerical result Fig. 

m 

The degree of mode coupling depends on the relative 
quantity between {Ei ~ EoY and coupling strength •q^VP' 
located in the second term of right-hand side of Eq. ([2|) . 
According to the relative quantity, the crossing point of il 
are determined and the transition of the type of ARC can 
be occur. It is a well-known phenomenon [iflll2;l20|. In 
Fig. [3]obtained from the toy model, (a) is about io(25, 9) 
and (b) is about io(25,8). Under the our coupling con- 
stant matched to the numerical result (Fig. [Ij , the point 
of transition from weak ARC to strong ARC for m = 2b 
is shown between these two IRs. 



IV. DYNAMICS OF 

MORPHOLOGY-DEPENDENT RESONANCES 

BY OPENNESS 

The dynamics from the numerical result and the toy 
model is explained in terms that the coupling between 
cavity itself and environment becomes strong as open- 
ness increases, and then the IRs are exquisitely affected 
by low-leaky ORs in large opening regime. Hence, a care 
should be taken to the identification of MDRs according 
to the history of ARCs. On a resonance position plot in 
the large opening regime, various IRs can coexist such as 
once and twice swapping resonances after strong ARCs, 
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FIG. 4: (Color online) Resonance positions at SI = 0.667 for 
TE polarization. Inverted triangles are 0(771, [ '"^•'" ] ) and the 
others are i(m,l) decided by considering ARC. Circles are 
\[m,l) changed from oo(m, [ "^^"^ ]). Dashed, dotted, and dot- 
dashed lines in (a) are Re[fc] = mO,, Re[fc] = m, and fc_B, 
respectively. The inset in (a) is for Q. — 0.5. Blue dotted line 
in (b) is 7. 



those superposed near weak ARCs and strong ARCs, and 
the uncoupled ones. Figure U] shows the resonance posi- 
tions in consideration of the ARCs. We included only 
the IRs and the low-leaky ORs for the sake of conve- 
nience. The dotted and the dashed lines in Fig. HJa) are 
Re[fc] = m and Re[fc] = mfi, respectively. In the effective 
potential analogy. Re [A:] = m corresponds to the top of 



the effective potential barrier, that is. the critical incident 
angle for IRs and Re[fc] = mil corresponds to the bottom 
of the barrier, that is, the incident angle 9 = tt/2. Effec- 
tive potential analogy is successfully applied to schematic 
description of the properties for IRs. For this reason, it is 
often used on the research about the resonances in dielec- 
tric disk. We, however, note that the effective potential 
analogy has some insufficient parts yet on applying to 
0Rs[3|. For a given m in Fig. [H^a), a circle is i(m, /) 
changed from Oo(to, [^^^^-^-j ) by the first strong ARC and 
an inverted triangle is o(m, [^^^^] ) changed from io{m, I) 
by the final strong ARC. The triangles between the cir- 
cle and the inverted triangle are i{m, I — 1) changed from 
io(m, I) by two consecutive strong ARCs. If we are able 
to get anyhow a correct position of o(m, [^^^^] ), it is nat- 
ural that io(TO, Z), whose Re[fcr]s are positioned between 
Oo(m, [^]) and o(m, [^]), are i{mj - 1) 0, be- 
cause 0(777., [^^^^^-'-j ) must be excluded from the ordering 
of I for 1(771, 1). 

Here, we newly use a line number for ORs defined as 
q = [22i±i] - /, (g = 0, 1, 2, • • • ) and, then, the low-leaky 
ORs are on the line of g = 0. The individual line of q 
except zero is generally parallel to the line of Re[/s] = 
771. However, the line of g = has a different slope. It 
surprisingly corresponds to the Brewster angle in large 
opening regime. In the viewpoint of semiclassics, the 
incident angle for an IR can be estimated by sin0 = 
777S1/Re[fc] and the Brewster angle is represented as 6b = 
arctan(Q). Consequently, one can obtain Re[fc] relevant 
to the Brewster angle, that is, ks = my/l + Q'^. In Fig. 
m^a) and inset, Re[fcr]s of low-leaky ORs are well-aligned 
on the ks line, whether they have experienced strong 
ARC or not. Also, by strong ARC in the opening regime 
where IRs and low-leaky ORs are coexistent in the same 
region of k plane, Im[A:r]s of low- leaky ORs settle down 
in the bottom region of the IR group as shown in Fig 
Hb). 

Figure [5lja) is the three-dimensional representation of 
the dynamics of io(25,/ < 13) and Oo(25,13) in Re[A:r], 
Ini[/sr] and Q space. For more detailed depiction and 
notation of reference marks, we projectively plot the dy- 
namics of Re[/ov] and Im[fcr] in Fig- ISH^) and Fig. Eljc), 
respectively. The dynamics of Re[A:r] for IRs and ORs 
in TM polarization exhibits almost linear and constant 
behaviors, respectively 0^. On the other hand, in TE 
case, the low-leaky ORs breaks the monotonic behav- 
ior of IRs as shown in Fig. \51[h) through the ARCs. 
As ri increases, Oo(25,13) interacts with io(25,/), where 
13 > I > 9, sequentially through a series of weak ARCs 
before it interacts with io(25,8) through the first strong 
ARC at SI ~ 0.481, i.e.. there is a transition from weak 
ARC to strong ARC. In the large opening regime, the 
point of strong ARC satisfies ks for given 777 and S7, 
since Re[fcr]s of IRs increase while Re[fcr]s of low-leaky 
ORs tend to maintain near ks- After the first strong 
ARC, Re[fcr] of io(25,8) stays near ks because the res- 
onance turns to 0(25, 13) by the exchange of their iden- 
tities. The strong ARCs between o(25, 13) and io(25,/). 
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FIG. 5; (Color online) Dynamics of MDR with m = 25 by 
openness. Solid lines and a dashed line are io(25,/ < 13) and 
Oo(25, 13), respectively. Dotted line in (b) and (c) is fcs and 
7, respectively. 



Q than the IRs for the same m, |Im[fcr]| of i(25,Z) in- 
creases by the interaction with o(25, 13). Moreover, the 
increased |Im[A:r]| after ARC is always below the bound 
line 7 due to the aftermath of the ARC. Hence, at a fixed 
value of 51, the IRs positioned below the 7 line are the 
resonances influenced by the ARCs. This is the reason 
why the IRs for TE polarization have generally lower Q 
than the corresponding modes for TM case [J, [5.:, 7, 8, 2a| 
and their imaginary part falls below the 7 line in the 
ordinary experimental regime of fl as shown in Fig|3l^b). 

Also, the dynamics gives a clear explanation for the 
breaking of the decreasing order of Q-factor in accordance 
with the increasing order of I near the Brewster angle. In 
a local region of complex k space under experimentally 
feasible circumstances, the grouping of resonances by the 
order of Q makes it possible for us to predict the char- 
acteristics of the resonances. If the effects of ARC are 
practically zero, the IR with higher I has generally low Q 
according to the effective potential analogy since a reso- 
nance with higher energy is less trapped in the effective 
potential @, S H, l25| ■ In TE case, however, the order of 
Q of IRs is reversed from a certain value of I due to ARCs 
with ORs as shown in Fig. HJc). For instance, at fi = 0.5, 
the order of Q corresponds to the order of I in TM case 
as (5i(25,9) = 85.876 > Qi(25,io) = 83.253 > (3i(25,ii) = 
82.528. On the other hand, in TE case, the order is 
reversed as Qi(25,9) = 27.116 < Qi(25,io) = 33.120 < 
Qi(25,ii) = 38.377. As f2 increases, this phenomenon is 
extended to the resonances with lower / due to the in- 
crease of the coupling strength. 

Such phenomena, only occur for TE polarization case, 
indicate that the existence of Brewster angle strongly ac- 
tivate the interactions between IRs and ORs as Q in- 
creases. In the dielectric disk for TE polarization, the IRs 
are sensitively affected by Brewster angle for the inter- 
nal cavity boundary. So, the IRs whose Re[fcJis near ks 
have larger leakage (|Im[/i:]|) than the others [2g. i.e., the 
more Re[A:r] of IRs approaches to /cb, the more the leak- 
age of the IR increases. An IR whose leakage is increased 
near the Brewster angle by the change of openness, as 
CI increases, approaches to a low-leaky OR in complex 
k space. Then, the coupling with the low-leaky OR be- 
comes relatively strong and the interaction between them 
is more activated. Such ARC for each IR occurs sequen- 
tially at other points of 57. Moreover, since the coupling 
strength becomes strong as fl increases, the transition to 
strong ARC from weak ARC occurs based on an IR with 
specific value of / as shown in Fig. [Sj 



where / ^ 8, sequentially occur at Re[fcr] ~ ks in the 
region of 57 ^ 0.481 until 57 becomes unity, where all 
resonances merge into the continuum. 

The dynamics of Im[fcr-] shows more interesting be- 
havior. Im[fcr] for a resonance is directly connected 
to the lifetime (r = — l/2clm[fc]) and Q-factor (Q = 
—Re[k]/2Iin[k]). On the space of Im[/i;,.] vs 57 as shown 
in Fig. [SJc), IR and low-leaky OR approach each other 
near the point of ARC. Since the low- leaky OR has lower 



V. SUMMARY 

Summarizing, we have demonstrated that the MDR 
dynamics by openness in a dielectric circular microdisk 
for TE polarization shows ARCs over two types of reso- 
nances due to the coupling with environment. The low- 
leaky outer resonances interact pronouncedly with the 
inner resonances through the opening channel near the 



Brewster angle. As a result, in the opening regime rele- 
vant to experiments, the resonance classification between 
inner and outer resonances is more complicated and the 
change of leakage for inner resonances occurs. The re- 
verse order of Q-factors of inner resonances which is dis- 
crepant from the effective potential analogy can be un- 
derstood in terms of the ARCs between inner and outer 
resonances. This study will be useful as a basis to ana- 
lyze intricate resonance dynamics in deformed dielectric 
microcavities. 
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